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Abstract 

We review the quark-meson coupling model, in which the quark degrees of freedom 
are explicitly involved to describe the properties of not only nuclear matter but also 
finite nuclei. Then, we present the electric and magnetic form factors for the proton 
bound in specific orbits for several closed-shell nuclei. 
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1 Introduction 

Whether or not quark degrees of freedom play a significant role in nuclei is one of the 
central questions in nuclear physics. Recently, tremendous efforts have been devoted to the 
study of medium modifications of hadron properties in a nucleus (which are in some sense 
precursors of the QCD phase transition) Jl[. At present lattice simulations have mainly 
been performed for finite temperature (T), with zero chemical potential. Therefore, many 
authors have investigated hadron properties at finite nuclear densities (p#) using effective 
theories [1-6]. 

In this paper, we will introduce one of these effective theories, namely the quark- 
meson coupling (QMC) model ||, £|, |5|, and then report the medium modification of the 
electromagnetic form factors of the bound proton in finite nuclei. 

2 The quark-meson coupling model 



2.1 Effect of nucleon structure 

Let us suppose that a free nucleon (at the origin) consists of three light (u and d) quarks 
under a (Lorentz scalar) confinement potential, V c . Then, the Dirac equation for the 
quark field ip q is given by 

[ry • d - m q - V c (r)]ip q (r) = 0, (1) 

where m q is the bare quark mass. 

Next we consider how Eq.(|l|) is modified when the nucleon is bound in static, uniformly 
distributed (iso-symmetric) nuclear matter. In the QMC model f|, |5| it is assumed 
that each quark feels scalar V£ and vector V£ potentials, which are generated by the 
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surrounding nucleons, as well as the confinement potential. Since the typical distance 
between two nucleons around normal nuclear matter density (po — 0.15 fm -3 ) is surely 
larger than the typical size of the nucleon (the radius Rn is ~ 0.8 fm), the interaction 
(except for the short-range part) between the nucleons should be colour singlet, namely 
a meson-exchange potential. Therefore, this assumption seems appropriate when the 
baryon density ps is not high. If we use the mean-field approximation (MFA) for the 
meson fields, Eq.(|l]) may be rewritten as 

[i 7 • d - (m q - V s q ) - V c (r) - 7 oKfH(r) = 0. (2) 

The potentials generated by the medium are constants because the matter distributes 
uniformly. As the nucleon is static, the time-derivative operator in the Dirac equation 
can be replaced by the quark energy — ie g . By analogy with the procedure applied to the 
nucleon in QHD |J, if we introduce the effective quark mass by m* = m q — the Dirac 
equation (|2|) can be rewritten in the same form as that in free space with the mass m* 
and the energy e q — instead of m q and e q . In other words, the vector interaction has 
no effect on the nucleon structure except for an overall phase in the quark wave function, 
which gives a shift in the nucleon energy. This fact does not depend on how to choose 
the confinement potential V c . Then, the nucleon energy at rest in the medium is given by 
En = M N (yv) + 3V£, where the effective nucleon mass M N depends on only the scalar 
potential. 

Now we extend this idea to finite nuclei. The solution of the general problem of a 
composite, quantum particle moving in background scalar and vector fields that vary 
with position is extremely difficult. One has, however, a chance to solve the particular 
problem of interest to us, namely light quarks confined in a nucleon which is itself bound 
in a finite nucleus, only because the nucleon motion is relatively slow and the quarks 
highly relativistic. Thus the Born-Oppenheimer approximation is naturally suited to the 
problem. Our approach in Ref. [|J was to start with a classical nucleon and to allow its 
internal structure to adjust to minimise the energy of three quarks in the ground-state of 
a system under constant scalar and vector fields, with values equal to those at the centre 
of the nucleon. Having solved the problem using the meson fields at the centre of the 
nucleon, one can use perturbation theory to correct for the variation of the scalar and 
vector fields across the nucleon bag. In first order perturbation theory only the spatial 
components of the vector potential give a non-vanishing contribution. This extra term is 
a correction to the spin-orbit force ||. 

As shown in Refs.0, ||, the basic result in the QMC model is that, in the scalar 
and vector meson fields, the nucleon behaves essentially as a point-like particle with an 
effective mass M N , which depends on the position through only the scalar field, moving 
in a vector potential. 

Let us suppose that the scalar and vector potentials in Eq.(|]) are mediated by the a 
and u mesons, and introduce their mean-field values, which now depend on position r, by 
V^{f) = g^er(f) and V£(r) = g^u(r), respectively, where g% (g^) is the coupling constant 
of the quark-cr (u) meson. Furthermore, we shall add the isovector vector meson p and 
the Coulomb field A(r) to describe finite nuclei realistically 0, Then, the effective 
Lagrangian density for finite nuclei, involving the quark degrees of freedom in the nucleon 
and the (structureless) meson fields in MFA, would be given by [|, |J 

- r N e 

Cqmc-i = # 7 -<9-M^-^u; 7 o-#p^6 7 o--(l+T3 )Aj ]ip 
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- + m> 2 ] + \\^f + m> 2 ] + \[{Vbf + m 2 /} + \{VA)\ (3) 

where if) and b are respectively the nucleoli and the p fields. m a , m w and m p are respec- 
tively the (constant) masses of the a, u and p mesons. g u and g p are respectively the cu-N 
and p-N coupling constants, which are given by g u = 3g q and g p = g q (where g q is the 
quark-p coupling constant). We call this model the QMC-I model ||. If we define the 
field-dependent cr-N coupling constant g ff (a) by 

M* N (a(r)) = M N - g a (a{f))a(f), (4) 

where is the free nucleon mass, it is easy to compare with QHD M. The difference 
between QMC-I and QHD lies only in the coupling constant g a , which depends on the 
scalar field in QMC-I while it is constant in QHD. However, this difference leads to a lot 
of favorable results. 

Here we consider the nucleon mass in matter further. The nucleon mass is a function 
of the scalar field. Because the scalar field is small at low density the nucleon mass may 
be expanded in terms of a as 



M* N = M N +[^pL\ a + Ll^JL) a 2 + -... (5) 




Since the interaction Hamiltonian between the nucleon and the a field at the quark 
level is given by H int = —3g q J df i/j q aijj q , the derivative of with respect to a is 
—3g q J df if; q ip q (= —3g q SN(cr)), where we have defined the quark-scalar density in the 
nucleon Sn(ct), which is itself a function of the scalar field. Because of a negative value 
of the derivative, the nucleon mass decreases in matter at low density. 

Furthermore, we define the scalar-density ratio Sn{ct)/ Sn(0) to be Cjv(ct) and the cx-N 
coupling constant in free space to be g a (i.e., g a = g„{a — 0) = g a — 3g^Sjv(0)). Using 
these quantities, we find 

M* N = M N - g a a - ^g a C' N (0)a 2 + ■ ■ ■ . (6) 

In general, is a decreasing function because the quark in matter is more relativistic 
than in free space. Thus, C' N (0) takes a negative value. If the nucleon were structureless 
Cn would not depend on the scalar field. Therefore, only the first two terms in the RHS 
of Eq.(||) remain, which is exactly the same as the equation for the effective nucleon mass 
in QHD §. 



2.2 Effect of meson structure 

It is true that not only the nucleon but also the mesons are built of quarks and anti- 
quarks, and that the mesons may change their properties in a nuclear medium. In Ref . || , 
we have studied the structure effects of both the nucleon and the mesons on the properties 
of finite nuclei. (We call this model QMC-II.) For further detailes on QMC-II, see Ref.0 
and references in the next section. 
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2.3 Infinite nuclear matter in QMC-I 

From the Lagrangian density Eq.(|^), we can easily find the total energy per nucleon 
Etat/A and the (constant) mean- field values of u and p (which are respectively given by 
baryon number conservation and the difference in proton and neutron densities). The 
scalar mean-field is given by a self-consistency condition, namely (dE tot /da) = 0. 

Now we need a model for the structure of the hadrons involved to perform actual 
calculations. We use the MIT bag model in static, spherical cavity approximation. In the 
present model, the bag constant B and the parameter (which accounts for the sum 
of the cm. and gluon fluctuation corrections 0]) for the nucleon are fixed to reproduce 
the free nucleon mass (Mjv = 939 MeV) and the free bag radius Rn = 0.8 fm. In the 
following we choose m q = 5 MeV, and set m a = 550 MeV, = 783 MeV and m p = 770 
MeV. (Variations of the quark mass and Rn only lead to numerically small changes in 
the calculated results |§.) We then find that B 1 ' 4 = 170.0 MeV and z N = 3.295. 

The coupling constants (g% and g%) are fixed to fit the average binding energy ( — 15.7 
MeV) at the saturation density (po) for symmetric nuclear matter. Furthermore, the p- 
N coupling constant is used to reproduce the bulk symmetry energy, 35 MeV. We then 
find |§: gl/A-K = 5.40, gl/Aix = 5.31 and gf/4ir = 6.93. Note that the model gives the 
nuclear incompressibility about 280 MeV, and the variations of the nucleon bag radius 
5R%/Rn = —0.02, the lowest eigenvalue Sx* n /xn = —0.16 and the root-mean-square 
radius of the quark wave function Sr*/r q = +0.02 at saturation density [|]. 

3 Applications 

The idea of the QMC model was first proposed by Guichon || in 1988, and later it has 
been developed by Adelaide group. In particular, Saito and Thomas have applied this 
model to various phenomena in nuclear physics. We here list up those applications: 

• Nuclear structure functions and the EMC effect 0, 

• Equation of state (EoS) for nuclear and neutron matter ||, 

• Charge symmetry breaking in nuclear matter - the Nolen-Schiffer anomaly 0, 

• Variation of hadron masses and matter properties in dense nuclear matter ||, [l~0|j , 



• Super-allowed Fermi beta-decay - the unitarity problem of the CKM matrix ele- 
ment HI, 



Properties of finite nuclei 0, [|, [12 



Naturalness in the QMC model |13 
Hyper nuclei [in 



• Electromagnetic form factors of the bound nucleon Jl5|, [16 

• In-medium Kaon properties IT? 



• Meson- nucleus bound states 118 
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Furthermore, the original QMC model was improved by Jin and Jennings |19 
which the bag constant B is allowed to decrease in nuclear matter. This version is called 
the modified QMC (MQMC) model. The MQMC model has also used to calculate the 
properties of finite nuclei, and the relationship among QMC, MQMC and QHD is clari- 
fied [Bo(l. 



4 Form factors of the bound nucleon 

As an example of recent applications, we show the in-medium modifications of the form 
factors of the nucleon [Tj|. In QMC, the quark wave function, as well as the nucleon 
wave function (both are Dirac spinors), are determined once a solution to equations of 
motion are found self-consistently M. The electromagnetic form factors for a proton 
bound in a specific orbit a, in local density approximation, are then simply given by 

G%AQ 2 ) = J G EM (Q 2 , p B {f))Ppa(f) dr, (7) 
where Ge,m(Q 2 pBif)) is the density-dependent form factor of a "proton" immersed in 



nuclear matter with a local baryon density, psif) ( see Ref.pofl). Using the calculated 
nucleon shell model wave functions, the local baryon density and the local proton density 
in the specified orbit a are easily evaluated. 

The notable medium modifications of the quark wavefunction inside the bound nucleon 
include a reduction of its frequency and an enhancement of the lower component of the 



Dirac spinor. As in earlier work ||21|| , the corrections arising from recoil and center of mass 
motion for the bag are made using the Peierls-Thouless projection method, combined with 
Lorentz contraction of the internal quark wave function and with the perturbative pion 
cloud added afterwards. Additional, possible effects of off-shell form factors and meson 
exchange currents are ignored in the present, exploratory study. The resulting nucleon 
electromagnetic form factors agree with experiment quite well in free space [21], at least 
for momentum transfer less than 1 GeV 2 . 

In order to reduce the theoretical uncertainties, we prefer to show the ratios of the 
form factors with respect to corresponding free space values. In Fig. [I] we show the ratios 
of the electric and magnetic form factors for 16 0, which has one s-state (lsi/2) and two 
p-states (lj?3/2 and 1^1/2)- As expected, both the electric and magnetic rms radii become 
slightly larger and the magnetic moment of the proton is also larger than the free value. 
The momentum dependence of the form factors for the ls-orbit nucleon is more supressed 
than those of the lp-states. This is because the inner orbit in 16 experiences a larger 
average baryon density. The magnetic moment for the ls-orbit nucleon increases by about 
7%, but, in the lp-orbits, it is reduced by about 2 — 3% from the ls-orbit value. The 
difference between two lp-orbits is rather small. For comparison, we also show in Fig. [I] 



the corresponding ratio of form factors (curves with triangles) using MQMC [19, where 
the bag constant is allowed to decrease by 10% The effect of a possible reduction 
in B is quite large and severely reduces the electromagnetic form factors of the bound 
nucleon. 

From the experimental point of view, the ratio, Ge/Gm, can be derived directly from 
the ratio of transverse to longitudinal polarization of the outgoing proton, with minimal 
systematic errors. We find that Ge/Gm runs roughly from 0.41 at Q 2 = to 0.28 at 
Q 2 = 1 GeV 2 for a proton in the ls-orbit in 16 O. The ratio of Ge/Gm with respect to the 
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Figure 1: Ratio of in-medium to free space electric and magnetic form factors for the ls- 
and lp-orbit nucleons of 16 O. The curves with triangles represent the corresponding ratio 
calculated in MQMC with a 10% reduction of the bag constant at po- 



corresponding free space ratio is presented in Fig. |2|. The result for the ls-orbit in 16 O is 
close to that in 4 He and about 2 — 3% lower than that for the lp-orbits in 16 0. We also 
find that the effect of the reduction in B has a significant effect on this ratio of ratios, 
especially for larger Q 2 . 

For completeness, we have also calculated the orbital electric and magnetic form factors 
for heavy nuclei such as 40 Ca and 208 Pb |16fl. Because of the larger central baryon density 



of heavy nuclei, the proton electric and magnetic form factors in the inner orbits suffer 
much stronger medium modifications than those in light nuclei. That is to say, the 
Q 2 dependence is further suppressed, while the magnetic moments appear to be larger. 
Surprisingly, the nucleons in peripheral orbits (1g?5/2, 2si/2, and for 40 Ca and 2rf 3 / 2 , 
l/in/2, and 3si/ 2 for 208 Pb) still endure significant medium effects, and comparable to 
those in 4 He. 



5 Conclusion 

In summary, we have reviewed the quark-meson coupling model, in which the quark de- 
grees of freedom are explicitly involved to describe not only the nuclear matter but also 
finite nuclei. Then, we have shown the electric and magnetic form factors for the pro- 
ton bound in specific orbits for several closed-shell nuclei. Generally the electromagnetic 
rms radii and the magnetic moments of the bound proton are increased by the medium 
modifications. In view of current experimental developments, including the ability to pre- 
cisely measure electron-nucleus scattering polarization observables, it should be possible 
to detect differences between the form factors in different shell model orbits. The current 
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Figure 2: Ratio of electric and magnetic form factors in-medium, divided by the free space 
ratio. As in Fig. [I], the curves with triangles represent the corresponding values calculated 
in MQMC. 



and future experiments at TJNAF and Mainz therefore promise to provide vital informa- 
tion with which to guide and constrain dynamic microscopic models for finite nuclei, and 
perhaps unambiguously isolate a signature for the role of quarks. 

I am pleased to thank Tony Thomas, Kazuo Tsushima, Ding Lu and Tony Williams 
for valuable discussions and comments. This work was supported by the Japan Society 
for the Promotion of Science. 
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